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ABSTRACT 

We present Doppler imaging and a Balmer line analysis of the weak-line T Tauri 
star TWA 17. Spectra were taken in 2006 with the UCL Echelle Spectrograph on the 
Anglo- Australian Telescope. Using least-squares deconvolution to improve the effective 
signal-to-noise ratio we produced a Doppler map of the surface spot distribution. This 
shows similar features to maps of other rapidly rotating T Tauri stars, i.e. a polar spot 
with more spots extending out of it down to the equator. 

In addition to the photosphcric variability, the chromospheric variability was 
studied using the Balmer emission. The mean Ha profile has a narrow compo- 
nent consistent with rotational broadening and a broad component extending out 
to ±220 km s _1 . The variability in Ha suggests that the chromosphere has at least 
one slingshot prominence 3i?* above the surface. 

Key words: stars: pre-main-sequence - stars: late-type - stars: chromospheres - 
stars: spots - stars: individual: TWA 17 - stars: magnetic fields 



The intrinsic variability and faintness of late-type pre-main 
sequence (PMS) stars has made it difficult to gain a com- 
plete picture of their development in the crucial first few 
tens of million years of their lives. In particular, this is 
because their variability often means that only long-term 
monitoring can give a useful picture of their behaviour. The 
tracks that young stars follow onto the main sequence are 
therefore not well defined observationally, and different the- 
oretical PMS tracks can predict widely differing parameters. 
Comprehensive observations of stars at each stage of their 
evolution to the main-sequence are needed. The discovery in 
the last ten years of nearby co-moving groups of young stars 
has provided s everal populations of sta rs with similar ages 
and distances (IZuckerman fc Sondl2004 ) from which targets 
suitable for extensive study can be identified. 

These groups include populations of T Tauri stars 
(TTS). Stars of spectral types F - M are considered to 
be T Tauri stars after they become optically visible. They 
are still contracting, and stable hydrogen burning has not 
started. The radiated energy of the star at this time is due 
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to its gravitational contraction. After the T Tauri phase has 
ended it may still be tens of millions of years before the star 
reaches the zero-age main sequence (ZAMS). During this 
time the star is in the post-T Tauri stage. 

A key feature of young stars is their magnetic activity 
which manifests itself, e.g., in the presence of large spots cov- 
ering a substantial proportion of the photosphere. Starspots 
are formed when magnetic flux tubes emerge through the 
surface of the star, reducing convection to these points and 
causing a reduction in temperature and the formation of 
cool (and therefore dark) spots. In non-accreting T Tauri 
stars the presence of cool spots is the dominant cause of the 
variability and gives rise to lightcurves with clearly defined 
periods (equal to the rotational period) and amplitudes of 
up to a f ew tenths of a magnitude in the V band. Doppler 
imaging dVogt fe Penrodlll983l : |P,icelll996l ; IStrassmeiedl2002l ; 

20081 ) can be used to produce maps of the spots on 
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rapidly rotating stars by taking advantage of the distortions 
the spots cause to rotationally broadened absorption lines. 

Prominences - material supported by the magnetic 
field lines extending from the surface, are also a magnetic 
phenomenon. In young stars 'slings hot' prominences have 
been observed, e.g . in Speed y Mic dDunstone et al.l 120061: 



peen observed, e.g . m speed y M ic (Uunstonc ct ai. zuuo; 
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ll98Sl : iDonati et al.lll999h and TWA 6 l|Skellv et al.ll2008h . 

These are clouds of material supported by the magnetic 
field at distances often beyond the corotation radius. The 
Balmer lines, in particular Ha, can be used to study the 
chromosphere, and slingshot prominences show up as in- 
creased emission outside the projected rotational velocity 
and as rapidly moving absorption transients inside the pro- 
jected rotational velocity. 



1.1 TWA 17 

The TW Hydrae Association (TWA) is a nearby co-moving 
group with an estimated age of 10 7 yr. Several of its mem- 
be rs have Hipparcos distan ces around 50 pc from Earth. 
In lLawson fc Crausel (|2005l ) two subgroups were identified 
within the TWA; the members numbered 1-13 lie at an av- 
erage distance of ~ 55pc, and those numbered 14-19 proba- 
bly belong to the Lower Centaurus-Crux (LCC) subgroup of 
the Ophiuchus-Scorpius-Centau rus OB-sta r asso ciation. A 
similar conclusion was reached bv lMamaielJ (|2005l ) using the 
moving cluster method to estimate distances to individual 
TWA stars. 

In lSkellv et al] (|2008h (hereafter Paper I) we discussed 
the weak-line T Tauri star TWA 6. Doppler images were 
presented and used to demonstrate that there is no differen- 
tial rotation in this star. Using the Balmer line emission 
variability we showed evidence for the presence of large- 
scale prominences and active regions that were cospatial 
with the spots on the photosphere. In this paper we anal- 
yse another reported member of the same association - 
TWA 1 7. TWA 17 has a spect ral type of K5, a period of 
0.69 d ([Lawson fc Crausel 120051 ), and a vsvai of 45 km s _1 
jWebb et al.lll999l) . 

The observations and data analysis of TWA 17 are set 
out in Sections [2] and [3] The Doppler images are presented 
in Section 2] followed by a discussion of the Balmer line 
emission in Section [5] In Section [6] we present arguments 
that TWA 17 in fact lies in the LCC subgroup. 



2 OBSERVATIONS AND DATA REDUCTION 

The observations, carried out 2006 February 11-20, are de- 
scribed in detail in Paper I. A brief summary is given here. 

Spectra were taken of TWA 6 and TWA 17 using the 
University College London Echelle Spectrometer (UCLES) 
at the Anglo- Australian Telescope (AAT). UCLES was cen- 
tred at 5500A giving a wavelength coverage of 4400 — 7200A. 
A slit width of 1.0" was used, giving a spectral resolution 
of 45000. As well as the target stars, we observed template 
stars of a variety of spectral types, B stars to obtain the tel- 
luric/sky spectrum, flats, darks, biases and ThAr arc frames. 
Observations were made 2006 February 11-13 and 17-20, 
but the second set of nights was concentrated on TWA 6 in 
order to measure differential rotation in that star. The full 
observing log is given in Tables 1 and 2 of Paper I. During 
the first set of nights 27 usable spectra of TWA 17 were ob- 
tained, and two more were taken during the second set of 
nights. 



3 LEAST-SQUARES DECONVOLUTION 

All data were red uced using the sta rlink routine echomop 
(more details in iMills etaL 1 1 19971 , and Paper I). To im- 
prove the effective signal-to-noise ratio (SNR) of the data we 
used the t echnique of least-sq uares deconvolution (LSD, de- 
scribed in lDonati et al.iri997h . This technique treats a spec- 
trum as a convolution of a line profile and a line-list (effec- 
tively a series of delta-functions). We produced a high SNR 
line profile from each spectrum by deconvolving the line- list, 
obtained from the Vienna Atomic Line Database (VALD, 
iKupka et al1l2q00h, using the program SPDECON (described 
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2500 absorption lines were used in the deconvolution. 
The mean SNR of the input spectra was 18, and of the LSD 
profiles was 630, giving a gain of 35. The continuum nor- 
malisation was successful in most cases. This was evident 
from the fact that in each profile the continuum regions at 
either side of the deconvolved profiles were at the same level. 
The exceptions were the two profiles produced from the ob- 
servations that took place during the second set of nights, 
indicating a small change in the instrumental set-up in the 
mean-time. This was rectified by fixing a slope to the lengths 
of continuum on either side of the profile and subtracting. 

In almost all of our spectra there was a strong absorp- 
tion feature close to km s~ (see Fig. [TJ. As discussed in 
Paper I, this is due to the presence of scattered moonlight 
contaminating the spectra. As TWA 17 has a lower usini 
than TWA 6 the 'Moon' profile covers a larger portion of the 
profiles. It is difficult to disentangle the size of this absorp- 
tion feature from the changes to the profile due to spots so 
we us ed the following strategy, described in lMarsden et all 
(2005). The errors in the region of the spectrum affected by 
the Moon contamination were increased by a factor of ten, 
although any large number would be adequate for this pur- 
pose. The profiles with the modified errors were then used as 
initial inputs for the Doppler imaging code. The pixels where 
the errors were increased will have less influence on the fits 
than the rest of the profile. Hence, the difference between the 
fits and the original profiles (the 'difference profiles') will be 
dominated by the Moon contamination. This contamination 
(which is essentially due to sunlight) was modelled using an 
unspotted, non-rotating star. A spectrum of such a reference 
star was deconvolved (using LSD) to provide a 'reference 
profile'. This reference profile is then fitted to each of our 
difference profiles. 

Fitting the difference profiles involved multiplying the 
reference profile by a different scale factor, as the level of 
contamination varies depending on the weather and Moon 
position. To find the scale factor a we minimise the value 
of £;(yi — azi) 2 for each profile, where y is the difference 
profile, a is the scale factor and z is the reference profile. 

In Fig. [2] the reference profile is shown. This was cre- 
ated by deconvolving a template spectrum (the K3 star 
HD 34673) and setting the flux equal to zero at velocities 
where the Moon is not considered to have an effect. The flux 
on the vertical axis has been rescaled to give a depth of 1.0, 
and the continuum has been shifted to 0.0. The result of this 
is that the reference profile can simply be multiplied by the 
appropriate factor and subtracted from the LSD profiles. In 
Fig. [3] the reference profile has been scaled in this way and 
overplotted on two example LSD profiles. 
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Figure 1. Line profile at phase 0.736. The data points are at 
intervals of 4 km s — 1 and the magnitude of the errors is indi- 
cated by the small bar on the bottom right. The large feature at 
km s _1 due to the Moon is obvious. 
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Figure 3. Scaled reference profiles (dashed lines) plotted over 
two of the TWA 17 LSD profiles. The values of a are 0.017 for 
the top plot and 0.033 for the bottom plot. 



Figure 2. The reference profile (see text) that has been multi- 
plied by the appropriate value of a before being subtracted from 
the LSD profiles for TWA 17. 



The values of a found for each profile are shown in Fig. 
[4] Two values of a are significantly higher than the others. 
This can be explained by the observing conditions of the 
individual spectra. The spectrum at phase 0.39 was taken 
when it was cloudy, and the spectrum at phase 0.67 had 
poor seeing. Finally, each profile had the difference profile, 
multiplied by the appropriate value of a, subtracted. The 
errors used at the velocities where the contamination was 
removed were then set equal to the size of the difference 
profile. 

It was noted that there is a radial velocity shift between 
profiles from different nights. This reached a maximum of 
7.6 km s" 1 between the first and last nights. This shift is due 
to changes in the instrumental set-up from night to night. 
We took account of this by cross-correlating between the 
telluric lines in different spectra and shifting the profiles by 
an appropriate velocity. 
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Figure 4. The scale factor a, plotted against phase. Each pro- 
file j has ajZ subtracted prior to reconstruction, where z is the 
'reference profile' plotted in Fig. [2] 
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Figure 5. Line profiles (1— sigma error bars) and fits (lines) for 
TWA 17. Errors have been increased where Moon contamination 
was removed. Phase and number of whole rotational periods since 
the first observation are shown to the right of each profile. 



4 DOPPLER IMAGING 
4.1 Stellar Parameters 

Of the 29 LSD profiles 28 were used as input for the 
Doppler imaging code DoTS (Doppler Tomography of Stars, 
ICollier CameronhoOll ). Details of DoTS are given in Paper I. 
One spectrum was discarded as it was the last one taken on 
the 2006 February 12 and was therefore noisy and badly 
contaminated by moonlight. 

First the stellar parameters, i.e., the wsini, period, in- 
clination angle and radial velocity were found. This was done 
by using DoTS to calculate x f° r a range of values; the ac- 
cepted value is where the minimum \ 2 ls found. Testing the 
code on artificial data has shown that this usually works 
even when the othe r parameters have not yet been found 
jBarnes et all [2000). However, the usini of the star must 
be determined at the same time as the equivalent width of 
the input spectral line. These two parameters were therefore 
found by simultaneously varying the v sin i and equivalent 




Figure 6. Greyscale images of the LSD profiles (left) and the 
corresponding Doppler imaging fits (right) of TWA 17, both di- 
vided by the mean profile. Spots appear as emission (white). The 
greyscale is given at the bottom. 



Table 1. 

TWA 17 



Previously-known 
References: [1]: 



and updated 



parameters for 
, , IZuckerman fc Sond (l2004h [2]: 
Cohen fc Kuhil J19791), [ 31 : iReidl ll2003Tl , [41 : lMamaiekl (I2005TI, [51: 
Laws on &: Craus el j2005l) . [6]: IWebb etalT l|l999h . [7]: ISiess et al] 



Parameter 


Previous 
Value 


New Value 


References 


Spectral type 


K5 


K5 


[1] 


Temp. [K] 




4400 ± 200 


[2] 


V magnitude 


12.7 




[3] 


V - R 


1.01 




[3] 


Distance [pc] 


90-163pc 


120 ± 40pc 


[1,4,5], this work 


usini [kms -1 ] 


45 ± 10 


47 ±1 


[3], this work 


^rad [ km s" 1 ] 


4.6 ±6 


11.1 ±0.1 


[6], this work 


Period [day] 


0.69 ± 0.01 


0.685 ± 0.001 


[5], this work 


Inclination [°] 




30 ± 10 


this work 


Radius [R©] 




i nc+0.16 
i - ZO -0.13 


this work 


Mass [M©] 




1.0 ±0.2 


[7], this work 


Lum. [L©] 




U ' M -0.19 


[3], this work 


Age [Myr] 


17 


12+ 9 
lz -6 


[5,7], this work 



width of the line, and finding the combination of values that 
gave the lowest v 2 - 

The inclination angle i is the most difficult parameter 
to constrain as large changes in i only yield small changes 
in x 2 - A value of 30° gave the lowest x 2 value with an error 
of ~ 10° on this value. Using the period, wsini and incli- 
nation angle gives a value of 1.251q \®Rq for the radius. 
The star can then be placed on pre-main sequence tracks in 
a Hertzsprung-Russell diagram of luminosity and tempera- 
ture. With an effective temperature of 4400K and a lumi- 
nosity of 0.531° '^ qL© the best-fit mass from the tracks from 



ISiess et al] l|2000D is 1.0±0.2M©. All resulting and adopted 
parameters are listed in Table. [1] some of which are taken 
from earlier work, in which case the reference is given. 

With these parameters it was then possible to fit the 
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data and produce Doppler images. Each pixel of the image 
is allowed to have a spot-filling fraction < / ^ 1, and 
the values of the / are allowed to vary as the code searches 
for the solution. For the photospheric temperature T p hot the 
pixels have / — > (/ 7^ to avoid singularities) . Pixels at the 
spot temperature T spo t have / = 1. During LSD the errors 
for the profiles are calculated using photon noise statistics, 
but this process can lead to und erestimation of the errors in 
the c ase of unpolarised profiles l|Donati fc Collier Cameron! 
1 19971 ). This leads to x 2 values greater than unity. A min- 
imum x 2 °f around 1.7 could be achieved but in this case 
we could not recover a maximum entropy image, so we re- 
laxed the final \ 2 requirement to 2.0. In Fig. [5]the data are 
shown (with the Moon spectrum subtracted) and the fits 
are overplotted. The increased error bars where the Moon 
spectrum has been subtracted can be seen clearly and it ap- 
pears that the fits are not affected by the Moon, as there 
are no spurious absorption features in the relevant velocity 
range. In Fig. [6] the data and fits are shown in greyscale, 
with a mean profile divided out in each case, and stacked in 
order of phase. The spot groups (the bright regions in the 
greyscale) are obvious. 

4.2 Surface Images 

The map resulting from the fits to the profiles (Figs [5] and 
[6]) is shown in Figs. [7] and [8] In common with other rapidly 
rotating young stars (including TWA 6) there is a polar 
spot and several spots at lower latitudes. The minimum and 
maximum temperatures indicated by the map are 3400K 
and 4400K respectively. This maximum temperature is de- 
termined using the spectral type (found by comparing with 
spectral standards). The minimum temperature is deter- 
mined by decreasing the temperature until the maximum 
filling factor falls below 1.0. 

The low SNR made it necessary to check that none of 
the features on the map are artefacts, caused by noise in 
the profiles. To do this the minimisation was repeated with 
random profiles removed. We wished to verify that the re- 
sulting image was similar even if we removed several profiles. 
However if too many profiles are removed it will leave large 
phase gaps. We elected to remove six randomly chosen pro- 
files each time. 

Once the profiles had been removed, we carried out the 
minimisation with 60 iterations and a target \ 2 °f 2-2. The 
target \ 2 na d to be slightly higher than the value of 2.0 
achieved earlier to ensure the we could reach the minimum 
X 2 each time. The procedure was repeated 30 times to pro- 
duce 30 separate images. The mean of the 30 images is shown 
in Fig. [9] with the normalised variance image overplotted as 
the contours. The contours line up well with the spots. This 
suggests that differences between the 30 images are due to 
differences in the exact intensity and shape of the spots in 
each image, rather than their locations. The only likely arte- 
fact is the feature at ~ 330° longitude and ~ 20° latitude 
which appears on a small scale in a number of images and 
thus also in the variance image. It appears in the map pro- 
duced using all data (Fig. [7} but is suppressed in the mean 
image. This feature is at a point where there is a gap in the 
phase coverage and we conclude that it is not real. 

The spot distribution is qualitatively similar to other 
maps of rapidly rotating PMS stars, i.e. there is a large 
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Figure 8. Polar projection of the spot map in Fig. [7] The dotted 
circles represent lines of latitude separated by 30° and the solid 
line shows the equator. Longitudes (corresponding to the longi- 
tudes in Fig. [7)1 are indicated by the numbers around the outside. 
Observation phases are indicated by tick marks. The greyscale is 
the same as in Fig. [jj 



polar spot and a lower spot coverage at all latitudes. The 
distribution can be c ompared with the flux tube models in 
iGranzer etafl (|2000h . The rotational velocity of TWA 17 is 
~ 50 times solar, and its mass is 1.0M®. The prediction for 
the spot map depends on the age we use for the star - 10 Myr 
for TWA or 17 Myr for LCC. In the former case TWA 17 
falls i nto the 'T' (T Tauri) classification in IGranzer et al.l 
(2000), and in the latter case the 'P' (pre-main sequence) 
classification applies. For the T case the models predict most 
of the emergence to be at the poles, with a smaller but non- 
zero spot coverage at lower latitudes. In the P case the spots 
are predicted to be in a band between 30 — 60° . 

Our recovered map looks more like the T case, where 
the star is assumed to have a very smal l radia tive core, con- 
sistent with the models of ISiess et all (|2000l) . These mod- 
els indicate that TWA 17 should have a small radiative 
core. However, previous Doppler images of post-T Tauri 
stars have also s hown high latitude and polar spots, e.g. 
RX J1508.64423 (|Donati et alj|2000l ) so as yet it is difficult 
to draw firm conclusions about the age of TWA 17 based on 
the spot distributions alone. 



4.3 Long-term evolution 

There is no contemporaneous photometry of TWA 17 but we 
can use our Doppler images to estimate the flux changes in 
the star as it rot ates and compare with ea rlier photometric 
measurements in lLawson fc Crausd l|2005h , whose observa- 
tions were carried out in May 2000. This is done by inte- 
grating the intensities across the visible disc of the star at 
each phase. The result is shown in Fig. 1101 The maximum 
change in brightn e ss of the star is just over 0.06 mag. In 
lLawson fc Crausel (|2005T ) a peak-to-peak AV of 0.12 mag 
was measured. Our lower estimate is reasonable given that 
Doppler imaging tends to underestimate the flux changes 
(we discuss this in more detail in Paper 1). The lightcurve 
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Figure 7. Mercator projection of the spot map produced using all the spectral data. Observation phases are indicated by the ticks at 
the bottom. The greyscale (given on the right) represents the filling factor. The mean spot filling factor at each latitude is also given. 
Phase increases from right to left, i.e., longitude of 360° is at phase and 0° is phase 1. 
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Figure 9. As Fig. [7]but showing the mean of all 30 spot maps (see text) in greyscale. The variance image is overplotted in red contours. 



in lLawson fc Crausel (|2005l ) is more symmetrical than Fig. 
[TOl indicating changes of the spot distribution in the inter- 
vening time. 

In Paper I we found that it was not possible to recover 
the A V of 0.49 mag that was measured in lLawson fc Crausel 
(|2005h . instead a value of 0.1 was calculated based on the 
Doppler images. We attributed this to a combination of a 
change in the spot coverage in the intervening years and a 
tendency of Dopple r imaging to underestimate flux changes 
|Unruh et al.ll 19951 ) . We do not find such a large discrepancy 
between the two TWA 17 values. This indicates that the 
fa ctor of 4-5 difference in the values for TWA 6 measured 
in lLawson fc Crausel (|2005l ) and our calculations cannot be 



accounted for by Doppler imaging underestimations alone. 



5 BALMER LINE ANALYSIS 

Balmer line emission is a tracer of the circumstellar envi- 
ronment and chromosphere of young stars such as TWA 17. 
The variability in Ha and the other Balmer lines will there- 
fore extend to velocities greater than vsini if the material 
is in corotation with the star. By studying the variation of 
these spectral lines it should be possible to identify events, 
such as flaring and prominences, in the chromosphere and 
corona. 



Doppler Images and Chromo spheric Variability of TWA 17 7 




0.4 0.6 
phose 

Figure 10. The V-band lightcurve predicted using the intensity 
map in Fig. [7] The symbols show the phases where spectra were 
taken. 
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5.1 Mean Balmer line profiles 

Figs [TT] and [12] show the mean Ha and H/3 profiles respec- 
tively. The mean has been taken of the continuum nor- 
malised profiles from each spectrum (also averaging the two 
orders that Ha and H/3 each appear in). A broadened ab- 
sorption spectrum has been subtracted off. The absorption 
spectrum was created by broadening the normalised tem- 
plate K7 spectrum to a «sini of 47 km s _1 (the usinz of 
TWA 17 found in Section gTTJ. Also plotted in Fig.[TT]is the 
normalised variance profile for the Ha profiles, given by 



V x = 



(n-l) 



1/2 



(1) 



where I^.i is the intensity of profile i at wavelength A, I\ 
is the mean of the intensities at wavelength A and n is the 
number of profiles. Equation [T] is based on th e expr ession 
for the variance profile given in Ijohns fc Basril |l995). The 
subtracted mean Ha profile has an equivalent width of 4.8A 
and a broad and narrow component. We consider the nar- 
row component to originate in the chromosphere. The origin 
of the higher velocities in the broad component will be dis- 
cussed later in this section. In Fig. [T2] both the subtracted 
and unsubtracted lines are shown because the subtraction 
makes a substantial difference to the shape of the line. The 
H/3 profile (with the broadened absorption spectrum sub- 
tracted off) has an equivalent width of l.lA. It can be fit- 
ted by a Gaussian profile with a FWHM of 70 km s^ 1 (also 
shown in Fig. I12[) . This is consistent with rotational broaden- 
ing and a v sini of 47 km s _1 us ing th e relat ionship between 
the FWHM and vsini given in lGravl (ll992T). 

Fo l lowing the example of, e.g.. |Petrov et all (|l994f ) and 
iHatzesl (1 19951 ) we attempted to fit the mean subtracted Ha 
profiles with a broad and narrow Gaussian. To do this we 
used x 2 minimisation to find the best fitting parameters (i.e., 
height, width and central position) for each Gaussian. The 
FWHM of the narrow and broad components found using 
this method are 75 km s _1 and 263 km s _1 respectively, 
giving a reduced \ 2 of 1.1. The fit is plotted in Fig. [13] The 
largest discrepancies between the fit and the data are at the 
positions of the telluric lines. The width of the narrow com- 
ponent is in agreement with the v sin i of 47 km s" 1 . There is 
a radial velocity shift of just under 8 km s _1 between the two 
components, which is slightly less than one pixel in resolved 
velocity (after binning), and we do not judge it significant. 
All parameters are listed in Table [2] 



Figure 11. The mean Ha line profile for TWA 17 is shown as the 
solid line. A rotationally broadened absorption profile (calculated 
as described in the text) has been subtracted. The normalised 
variance profile for the Ha line is shown as the dashed line. The 
vertical axis on the left-hand side refers to the mean Ha and the 
right-hand axis refers to the variance profile. The small dips in 
both profiles are due to telluric lines. 
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Figure 12. Mean H/3 profile of TWA 17 before (top) and af- 
ter (bottom) a rotationally broadened K5 template spectrum has 
been subtracted off. The broadened template spectrum is shown 
as the dot— dashed line. The dashed line shows the Gaussian fit to 
the subtracted profile. 



5.2 Variance profile 

The variance profile has a full-width at half-maximum 
(FWHM) of 178 km s _1 and indicates that there is vari- 
ability out to around ±260 km s~ . There a slightly higher 
variability on the blue side of the emission. However, this 
asymmetry is not apparent in the mean Ha profile. The 
narrow component of the variance profile is wider than the 
narrow component of the mean profile. As the wings may 
vary more than the line centre the full-width at zero in- 
tensity (FWZI) may be the more relevant parameter. We 
estimate that the FWZI is 260-280 km s" 1 although it is 
difficult to determine due to the broad 'pedestal' compo- 
nent. The normalised variance profile for the H/3 profile has 
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Figure 13. The solid line shows the mean Ha profile. The dashed- 
dotted lines are the two Gaussians that were fitted to the mean 
profile. They have FWHM of 75 km s" 1 and 263 km s" 1 . The 
dashed line is the sum of the Gaussians. This fit to the line has a 
reduced \ 2 of 1.1. 



Table 2. Parameters of Gaussians used to fit Ha lines, for the 
narrow and broad components. 

Component Amplitude Centroid Width 

[ km s -1 ] [ km s" 1 ] 
Narrow L3 LO 75 
Broad L0 -7J3 263 

essentially the same shape as the mean profile but is quite 

noisy and is not shown. 

In previous work, e.g., iFernandez et alj l|2004 ). the 
broad 'pedestal' part of the emission has been explained by 
microflaring, i.e., small flares that cannot be resolved indi- 
vidually. The motions of the material within these flares can 
explain the emission seen in the wings of the line profiles. 
Microflaring may therefore be a contributor to the variabil- 
ity in the Ha line. 

5.3 Variability in the Balmer lines 

In Fig. [I4]the greyscale image of the Ha variability is shown. 
This was produced by stacking all the Ha emission lines in 
order of phase and subtracting the mean Ha profile (Fig. 

The greyscale at each pixel therefore represents the dif- 
ference in flux between that pixel and the mean, hence a 
greyscale value of zero indicates the emission at that pixel 
is equal to the mean. The velocity scale in Fig. [14] has been 
shifted to take account of the radial velocity. The asymme- 
try that was evident in the normalised variance profile (see 
Fig. is also visible in Fig. [14] as the left (blue) side of the 
image has more features (i.e bright or dark regions) than the 
right (red) side. The equivalent H/3 image was very noisy (as 
it is at the blue end of the spectrum where there are fewer 
counts) and therefore did not show any clear features. Con- 
sequently the H/3 image is not shown here. 

On TWA 6 we observed signatures of active regions that 
appeared to be cospatial with the spots found in the Doppler 
imaging. However, there is no evidence of such an alignment 
between the spots and active Ha regions on TWA 17. 
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Figure 14. Ha profiles of TWA 17 stacked, with a mean profile 
subtracted, shown in greyscale. The scale bar is shown at the 
bottom. The vertical lines indicate ±-usini. The dash-dot curve 
is a sinusoid indicating the position of the feature identified as a 
slingshot prominence (see text). The vertical axis shows phases 
from - 1.05, in order to show the continuation of the bright 
feature. 



Variability in Ha can be due to scattering of Ha photons 
from the stellar surface by circumstellar hydrogen clouds. If 
these clouds are in front of the stellar disc, photons will 
be scattered out of the line-of-sight. This leads to reduced 
emission, which will appear in images such as Fig.[l4]as dark 
regions. On the other hand if the clouds lie off the stellar 
disc photons will be scattered into the line-of-sight leading 
to increased emission. Variability at velocities above i;sini 
can either be caused by high-lying material in co-rotation 
with the star (e.g., slingshot prominences, see Section [TJ, 
or fast-moving material closer to the surface, e.g. related to 
flares. 

If a slingshot prominence is present and observed at all 
phases it should be possible to fit a sinusoid to the signa- 
tures. It may pass in front of the stellar disc, depending on 
its height and inclination, in which case it will most likely 
cause a decrease in emission. This will appear as a dark di- 
agonal band going from blue to red as the material is in 
corotation with the star. Scattering of Ha photons can also 
be caused by flare material, and infalling gas, which may 
also be accompanied by velocity changes that can indicate 
the direction of motion of the material. 

Studying Fig. [14] we see two features that may be ex- 
plained by these phenomena. The first of these is a bright 
region between ±80 km s _1 and phases ~ 0.4 — 0.8 going 
from red to blue, which may indicate that it is material on 
the far side of the star that does not go out of sight due 
to the inclination of the system. The second feature is a 
prominent brightening between 0.65 and 1.075 that extends 
between -220 km s -1 and approximately +120 km s _1 . It 
is possible that this is a slingshot prominence of the type 
discussed above, although this is difficult to confirm as we 
do not observe it for a whole rotation period. In Fig. 1141 this 
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brightening is overlaid with a sinusoid that represents the 
expected path of a prominence of this size. As the obser- 
vations that gave these spectra were all taken in one night 
(2006 February 13) we have no handle on its lifetime. It does 
not appear in absorption between iusini, which may indi- 
cate that it is indeed a prominence that lies at a distance of 
more 3i?» above the surface. 

In summary, although microflaring may be a contribu- 
tor to the variance, the large scale features that have high 
velocities in Fig.[l4]are important contributors to the overall 
variability in TWA 17. 



6 DISTANCE TO TWA 17 AND THE 
TW HYDRAE ASSOCIATION 

In this section we compare our findings for TWA 6 and 17 
and consider the question raised in Section[T]of whether they 
have similar distances and ages and hence can be considered 
members of the same association. 

Using the parameters found in Section 14.11 the radius 
of TWA 17 was estimated as I^S^q^Rq compared with 

1.05io;igR Q for TWA 6 - If TWA 6 and TWA 17 were to 
belong to the same association and were thus at a similar 
distance, we would not expect TWA 17 to have a higher 
effective temperature and be fainter than TWA 6. 

However, this discrepancy can be resolved if we consider 
that TWA 17 may be more distant than TWA 6. For the pa- 
rameters determined in this paper (specifically the radius) 
to be consistent with the spectral type and luminosity mea- 
sured previously a distance of 190±130pc would be required, 
where the error is dominated by the error in the radius that 
has been propagated. Given that TW Hya has a Hipparcos 
distance of 56 pc this might suggest that TWA 17 lies in the 
LCC subgroup, although the magnitude of the errors makes 
it difficult to draw firm conclusions on this point, 

Previous estimates for TWA 17 have suggested 
distances between 100 and 150 pc. For instance, in 
lLawson fc Crausel (2005) two separate populations have 
been identified within the stars previously identified as mem- 
bers of the TWA, those numbered 1-14 are the 'true' 
TW Hydrae association and the others are members of the 
more distant Lower Centaurus-Crux (LCC) subgroup of the 
Ophiuchus-Scorpius-Centaurus association, at an average 
distance of 90 p c. An even more distant estimate is given in 
iMamaiekl (120051) where the distanc e to T WA 17 is 163 ± 46 
pc. In addition [Zuckerman fc Sond l|2004h give a luminosity 
distance of 133 pc for TWA 17 and conclude based on this 
that it is a member of the LCC. The results from this paper 
are consistent with this interpretation. 

Using the luminosity and temperature values which we 
have adopted for these two stars we have placed them in a 
Hertzsprung-Russell (HR) Diagram in Fig. 1151 along with 
the brightest star in the association: TW Hya. The parame- 
ters for TW Hya come from p revious work using a spectral 
type of K8 (iTorres et al.ll2000l ). and a lu minosity of 0.24L fl 
(|Yanget al.l 120051 ). The tracks are from ISiess et all (|2000h . 
All three stars fall closely below the 10 Myr isochrone. 

We have looked at several other age indicators for the 
two stars. Both have Li 6707A absorption. TWA 6 has a 
Li 6707 A equivalent width of 0.6 A and TWA 17 has an 
equivalent width of 0.5A, suggesting that TWA 6 may be 
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Figure 15. Luminosity against effective temperature for TWA 6, 
TWA 17 and the brightest star in the association, TW Hya. The 
solid lines are model evolutionary tracks in units of solar masses. 
The dashed lines show isoch rones annotated in units of millions 
of years. The tracks are from lSiess et aL (2000). 

slightly younger. Curiously, we measure a Li EW of 0.35A 
for TW Hya from our reliable spectra which would imply 
that TW Hya is older. However, the lower equivalent width 
in TW Hya might be a consequence of photospheric veiling 
due to disk accretion. The difference in effective temperature 
can also lead to a spread in the lit hium EW of stars of 
identical ages l|Soderblom et al.lll993D . 

TWA 6 and TWA 17 both have Ha equivalent widths 
of around ~ 5A which classifies t hem a s WTTS according to 
iBarrado v Navascues fc Martini l|2003l ) . The Doppler images 
of both stars are similar, both have a polar spot and spots 
at all latitudes to the equator. As stated in Section 14.21 this 
corres ponds most closely with the T Tauri models in lGranzerl 
( 2004 ) (rather than post-T Tauri or ZAMS), but it is not 
possible to draw firm conclusions on the star's age based on 
this. 

The LCC stars are believed to be older than the TWA. 
Although these results presented here appear to show that 
the stars are the same age we cannot rule out the possibility 
that TWA 17 is older than TWA 6 as the error bars in Fig. 
1151 encompass a wide range of ages. 



7 DIFFERENTIAL ROTATION 

Due to the observing conditions in the second half of the ob- 
serving run the observations concentrated on TWA 6, mak- 
ing it difficult to measure differential rotation on TWA 17. 
We can nevertheless use the Doppler imaging code with the 
differential rotation parameter AO as an input, and use \ 2 
minimisa tion to attempt to place an u pper limit on its value 
(see e.g.. ICollier Cameron et al1l2002l ). 

When reduced \ 2 i s plotted as a function of angu- 
lar velocity Q and Afi, the resulting surface is noisy. The 
minimum contours form an extended ridge that crosses the 
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Q — Afi surface diagonally and encompasses a large range of 
differential rotation values. It is thus difficult to determine 
the confidence limits for the differential rotation of TWA 17. 
This is not surprising as we have only two profiles in the lat- 
ter half of the run, and the phases do not correspond to any 
of the phases in the earlier part of the observing run. Fur- 
thermore, the relatively low inclination of the star (~ 30°) 
reduces the resolution at lower latitudes allowing a wider 
range of differential rotation values to fit the data. The pe- 
riod quoted in Table Q] assumes zero differential rotation. 

As stated in Section[U] TWA 17 should be an interesting 
candidate for a more accurate differential rotation measure- 
ment. This is due to its position on the HR diagram (Fig. ll5p 
at the turn-of f from the Hayash i to the Henyey tracks, ac- 
cording to the Sicss et all (|2000h models. 



8 CONCLUSIONS 

We have presented Doppler images of TWA 17. A compar- 
ison with TWA 6 (Paper I) shows broadly similar features, 
i.e. a polar spot and further spots extending out of it down 
to the equator. The latitudinal distributio n of the spots is 
similar to that predicted by the models in iGranzerl (|2004l ) 
for a rapidly rotating T Tauri star. 

Despite the presence of substantial absorption features 
in the spectrum due to contamination by Moonlight, and 
relatively low SNR we are confident that the map obtained is 
robust, as the minimisation was repeated a number of times 
with random phases removed, giving similar maps each time. 

We also studied the chromospheric emission. The mean 
H/3 emission is consistent with rotational broadening, once 
a template absorption spectrum has been subtracted off. 
Ha has a broad component in addition to the narrow com- 
ponent. The presence of a broad component may indicate 
micro-flaring, although this cannot explain all of the vari- 
ability in the Ha profiles. The Ha line shows variability 
at velocities greater than usini, indicating that there is 
fast-moving material in the chromosphere, possibly flares 
or prominences. There is at least one feature in the emission 
that may be caused by a prominence. It lies at a distance 
of over 4_R* from the axis of rotation. However it is only 
apparent in one night's observations, making it difficult to 
draw conclusions about the lifetime. In Paper I we found ev- 
idence for chromospheric active regions that are co-spatial 
with the photospheric spots on TWA 6, but this work does 
not find any evidence for such regions on TWA 17. 

We obtain a distance of 190 ± 130 pc for TWA 17, which 
suggests that it does not lie in the TW Hydrae association 
but in a more distant subgroup. However, the age estimators 
indicate that TWA 17 and the bona fide association member 
TWA 6 are approximately of the same age. 

In the Balmer line analysis, there are signatures of high- 
lying prominences with lifetimes of several rotational pe- 
riods. In Paper I at least one prominence was visible on 
TWA 6. The increased emission at velocities outside ±usini 
suggest that similar features exist in TWA 17. The fact that 
we see evidence of such features in both stars would justify 
a dedicated survey of a young stellar cluster to gain a better 
picture of the frequency and nature of slingshot prominences 
in T Tauri stars. 
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